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Abstract 

The growth rate of matter perturbations is computed in a magnetized environment for 
the ACDM and wCDM paradigms. It is argued that the baryons do not necessarily follow 
into the dark matter potential wells after they are released from the drag of the photons. 
The baryonic evolution equations inherit a forcing term whose explicit form depends on the 
plasma description and can be deduced, for instance, in the resistive magnetohydrodynam- 
ical approximation. After deriving an analytical expression for the growth rate applicable 
when dark energy does not cluster, the effects of relativistic corrections and of the inhomo- 
geneities associated with the other species of the plasma are taken into account numerically. 
The spectral amplitudes and slopes of the stochastic magnetic background are selected to 
avoid appreciable distortions in the measured temperature and polarization anisotropics of 
the Cosmic Microwave Background. The growth of structures in the current paradigms 
of structure formation represents a complementary probe of large-scale magnetism in the 
same way as the shape of the growth factor and the associated indices can be used, in the 
conventional lore, to discriminate between competing scenarios of dark energy or even to 
distinguish different models of gravity. 
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1 Formulation of the problem 



The growth rate of matter inhomogeneities is customarily employed to distinguish the physi- 
cal features of structure formation scenarios. In terms of the density contrast of the inhomo- 
geneities associated with pressureless matter (i.e. 5 m (x,r)) the growth rate will be defined 
as [1J: 

f(x,y) = j m , S m (x,r) = — 8 c (x,r) + — 8 h (x,r), (1.1) 
ay cjmo wmo 

where y will denote, conventionally, the natural logarithm of the scale factoi@, while 5 C and 

<5b are, respectively, the cold dark matter (CDM) density contrast and the baryon density 

contrast; furthermore, following the standard notations, Uxo = h^Vlxo where Qxo is the 

(present) critical fraction of the species X of the plasma. 

The growth rate depends on the dynamical features of the dark energy background and 

it is affected by relativistic corrections which become relevant for typical wavelengths of the 

order of the Hubble radius (2j |3] . A rather useful parametrization of the growth rate involves 

the so-called growth index denoted hereunder with 7: 

f( y ) = n&, n M = ^, (1.2) 

Pt 

where f2 M = Qm(v) is n °t simply evaluated at the present time but it is a y-dependent 
quantity. In spite of the homogeneous parametrization of Eq. ( 11. 2ft . the density contrast 
of Eq. ( 11.11) is in general inhomogeneous (see e.g. [21 13]) and its spatial dependence ranges 
from wavenumbers comparable (and possibly smaller) than the equality wavenumber fc cq to 
wavenumbers directly probed by the currently available large-scale structure data, i.e. from 
^min = 0.01 /i Mpc _1 to, approximately, fc max = 0.25 h Mpc -1 . The range /c min < k < fc max 
includes also the scale at which the spectrum becomes nonlinear, i.e. k n \ > 0.09 ho Mpc -1 . 
According to the WMAP 7yr data pQElE], k eq = 0. 0097 A°:^ \ Mpc -1 corresponding to an 
effective equality multipole £ eq = 137.5 ± 4.3. Typical wavenumbers k ^> k cq crossed inside 
the Hubble volume before matter-radiation equality. Conversely the very large length-scales 
(relevant for the region of the Sachs- Wolfe (SW) plateau) fall into the complementary regime 
k <C k cq . The SW contribution peaks for comoving wavenumbers of the order of k sw = 
0.0002 Mpc -1 while the integrated SW contribution (typical of the ACDM paradigm and of 
its extensions) peaks between 0.001 Mpc" 1 and 0.01 Mpc -1 (see, e.g. [7J El U QH HI H2] ) • 
The definitional of 7 adopted in Eq. ( II .2p (see, e.g. (131 HH [15]) implies, in the context 



2 The present value of the scale factor is normalized to 1, i.e. ao = 1; moreover we will conventionally 
denote y = In a where a = a/a<j c and Ode denotes the value of the scale factor at which the critical fractions 
of matter and dark energy are equal. To avoid potential confusion we stress that, in the present script, the 
notation "In " is employed for the natural logarithm while the symbol "log " denotes the common logarithm 
of the corresponding quantity. 

3 While different parametrizations of the dark-energy component may lead to slightly different semi- 
analytic determinations of the growth index, these differences will not be essential for the purposes of the 
present investigation. 
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of the vanilla ACDM scenario, that 7 ~ 0.55. If the dark-energy component is characterized 
by a constant barotropic ratio Wde = Pde/Pde the growth index shall depend on Wde- The 
reference model where Wde is constant is customarily dubbed wCDM (see, e.g. [U El E]) and 
this terminology will also be employed in the forthcoming sections of the present paper. 

It is not implausible that, in the next decade, dedicated observations will be able to 
probe the growth rate and, more generally, the growth of structures either in the ACDM 
scenario or in one of its neighboring extensions such as the wCDM scenario (see, e.g. [T6|[T7] 
for two dedicated reviews on models of dark energy). An observational scrutiny could even 
allow for direct cosmological tests of gravity theories as convincingly argued in [TBI [19] . At 
the moment the large-scale observables directly or indirectly related to the growth rate of 
matter perturbations are the integrated Sachs- Wolfe effect [7J El El QUI Ell EE2J probed by 
Cosmic Microwave Background (CMB) experiments and the large-scale galaxy distribution 
scrutinized, for instance, by the Sloan digital survey [201 Ell [22] • Useful complementary 
informations can be also obtained from x-ray bright clusters [231 121], from Lyman-a power 
spectra and from weak lensing observations [231 123] . 

Large-scale magnetic fields can affect the growth rate of matter perturbations. Conse- 
quently the determinations of the growth index can be used as a potential probe of large-scale 
magnetism. The aim of the present investigation is to initiate a more systematic scrutiny of 
the effects of large-scale magnetism on the growth of structure in the same way as, in recent 
years, the effects of large-scale magnetic fields on the different observables related to the tem- 
perature and polarization anisotropies of the CMB have been analyzed (see [121 123 1211 [29] 
and references therein). 

According to the WMAP 7yr data [U E] photons decouple at a typical redshift Zd ec = 
1088.2 ± 1.2. Shortly after, also baryons are freed from the Compton drag of the photons at 
a redshift Zd TSug = 1020.3 ± 1.4. Within a simple magnetohydrodynamical appproximation, 
for z < Zdrag the baryon fluid is still strongly coupled and globally neutral while the photons 
follow, approximately, null geodesies. Both for z > Zd ec and for z > Zd ra ,g the equations 
governing the evolution of the baryon density contrast inherit a direct dependence on the 
magnetic field. The total density contrast 5 m , being the weighted sum of the CDM and of the 
baryon density contrasts, will also be affected by the large-scale magnetic field. The growth 
rate appearing in Eq. (II. ip . the growth index of Eq. (11. 2p will therefore bear the mark of 
large-scale magnetism. As direct bounds on the parameters of the magnetized ACDM and 
wCDM scenarios have been obtained from the analysis of the temperature and polarization 
anisotropies [T21 [27], it is conceivable that complementary bounds can be deduced from 
the study of the growth of structures. It is nonetheless clear that the impact of large-scale 
magnetism on the growth of structures must me carefully computed beyond the conventional 
perspective where the evolution equation of 5\> is forced by the evolution of 5 C so that, 
asymptotically, 5^ will approach the value of 5 C : this phenomenon is dubbed, at least in the 
linear regime, by saying that the baryons fall into the dark matter potential wells. This is, in 
practice, also the situation occurring in the ACDM paradigm where the dark energy does not 
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cluster and the density contrast of the dark energy vanishes (i.e. 5d c = 0). Such a perspective 
is quantitatively and qualitatively challenged by the presence of large-scale magnetic fields. 
In the latter situation the evolution of the baryons is the result of the competition of the 
dark-matter forcing with the magnetohydrodynamical Lorentz force. 

The implications of large-scale magnetic fields for the paradigms of structure formation 
have a long history dating back to the contributions of Hoyle [3U] and Zeldovich [21] (see 
also [32J E3] for two dedicated reviews on large-scale magnetic fields in astrophysics and 
cosmology). Peebles (see, e.g. [HEl]) and Wasserman [35] argued, within a non-relativistic 
treatment, that the (comoving) coherence scale of the magnetic field should be larger than 
(but of the order of) the comoving magnetic Jeans length_| 



if the magnetic field is expected to cause an appreciable contribution to the total matter 
density contrast after photon decoupling. The point of view implicitly endorsed by [U 135] 
was, however, to seed structure formation in a purely baryonic Universe and without the 
help of any adiabatic or isocurvature initial conditions. The ideas spelled out in [U [35] 
have been partially revisited in a pure CDM scenario and in a non relativistic set-up by 
Coles [36] without dark energy. The perspective of the present paper is not to endorse 
an alternative structure formation paradigm but rather to scrutinize the impact of post- 
recombination magnetic fields on the current ACDM and u>CDM scenarios. In doing so it 
will be essential to profit of the results obtained from the analysis of the magnetized CMB 
anisotropies [12j |27]. To prevent the distortions induced by a stochastic magnetic field on 
the temperature and polarization anisotropies the parameters of the magnetic power spectra 
must be in specific ranges which slightly differ in the case of the ACDM and in the case of 
the wCDM scenario. In summary the analysis of post-recombination effects of the large-scale 
magnetic fields in the ACDM and wCDM paradigms is still an open problem as the lack of 
specific predictions for the magnetized growth rate shows; one of the purposes of the present 
paper is to remedy this situation by setting a general framework where these themes can be 
quantitatively addressed. 

The layout of the present paper can be spelled out as follows. In sec. [2] the govern- 
ing equations of the system will be introduced in the synchronous gauge; the connection 
of the discussion with other gauge-dependent and gauge-invariant treatments will be briefly 
outlined. In sec. [3] the computation of the growth factor shall be specifically addressed. In 
section H] the analytical estimates will be corroborated by a specific numerical analysis includ- 
ing the effect of the appropriate relativistic corrections and of the dark energy fluctuations 
on the final shape of the growth rate. The final remarks and the concluding perspectives are 
collected in section [51 

We use here the notations employed in the bulk of the paper and Bl denotes the comoving amplitude 
of the regularized magnetic field measured in nG. Recall that 1 nG = 10~ 4 nT = 6.9241 x 10~ 29 GeV 2 . 
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2 Governing equations 



2.1 ACDM and wCDM backgrounds 

The toy models extending and complementing the vanilla ACDM scenario are classified in 
terms of a set of reference paradigms relaxing at least one of the assumptions of the conven- 
tional ACDM model. For instance, while in the ACDM case the dark energy perturbations 
are totally absent, in the wCDM (where w stands for the barotropic index of the dark en- 
ergy background) the dark energy fluctuations affect, indirectly, the evolution of the matter 
density contrast. Both in the ACDM and in the u>CDM models the background geometry is 
described by a spatially flat metric of the type g^^r) = a 2 {r)rj^ (where rj^ is the Minkowski 
metric). The Friedmann-Lemaitre equations take their standard form 

n 2 = ^^p t , H 2 -n> = 4nGa 2 ( Pt + p t ), H = -, (2.1) 
3 a 

where the prime denotes a derivation with respect to the conformal time coordinate r. The 

explicit form of the total energy density and of the total enthalpy density appearing in Eq. 

(12. ip is given by: 

Pt = Pc + Pi + Pc + P 7 + Pv + Pde, (2.2) 

4 

Pt + Pt = Pe + Pi + Pc + ~{Pu + P 7 ) + (W dc + I) Pde] (2.3) 

the subscripts in Eqs. (12. 2R and (12.31) denote, respectively, the electrons, the ions, the CDM 
species, the photons, the neutrinos (assumed massless both in the ACDM and in the wCDM 
paradigms) and the dark energy contribution. For z < z dec the baryonic matter density 
Pb = Pe + Pi is customarily introduced where p e = m e n and p\ = m;n are, respectively, 
the electron and ion matter densities. The comoving concentrations of the electrons and of 
the ions (i.e. no = a%) coincide becanse of the electric neutrality of the plasma B In the 
ACDM (and wCDM) paradigms it is sometimes practical to define further combinations of 
the energy and matter densities: 

Pm = Pe + Pi + Pc = Pb + Pc, Pr = P^ + P 7 , (2.4) 
implying that Eqs. (12.11) and (I2.2p can also be written as^|: 



H 2 — "Hq ^ M0 
Z(a) = — - = 

K J n 2 22 



I _|_ [ ^£5.^ _|_ ( ^£ x 
a 



w de (1 - Q M ) + (- - w dc ^JQ R 



(2.5) 
(2.6) 



5 The plasma is globally neutral rii = n c = uq and the common value of the electron and ion 
concentrations can be expressed as no = rj-^n^ where n 7 is the comoving concentration of photons, 
?/ b = 6.219 x 10- 10 [w b o/(0.02773)][T 7 o/(2.725K)]- 3 and u hQ is, as in Eq. CCLl>> the critical fraction of 
baryonic matter multiplied by h^. 

6 Notice that since the scale factor is normalized in such a way that ag = I, then "Ho = Hq = 
100 h km/[secMpc]. 



5 



where the redshifts of matter radiation equality (i.e. a eq ) and dark energy dominance (i.e. 
ad e ) are defined in terms of the present critical fractions of matter, radiation and dark energy: 

z eq + 1 = — = — , z dc + 1 = — = — = I — 1 . (2.7) 

«eq "RO «de ViideO/ \"do / 

In Eq. (12. 7p fi M0 , ^ro an d ^deo denote, respectively, the present critical fractions of matter, 
radiation and dark energy. Conversely, in Eqs. (12.51) and ( 12.61) the time-dependent critical 
fractions are defined as: 



(,2 cq + l)a 3u, dc + ( Zdc + l)o;3^ de -l + 



cq 



feq + i: 



nde{a) (z eq + l)a^ + (z de + 1)^-1 + ( 2cq + 1) ' (2 - 8) 

where, as in Eq. (II. ip . a = a/a de and fl^a) = [1 — fi M (a) — f2 de (a)] which also implies 
af2 M (ct) = (zdc + l)^R(tt)/(^eq + l)- When a ^> a eq the contribution of ^(a) is by definition 
subleading in comparison with Qm(&) and fid e (a). In the limit £7r <C ^m, ^(«) can be 
approximated as 

lim Z(a) Z(a) = - J - ^ de (l - Q M ). (2.9) 

By introducing the natural logarithm of the normalized scale factor, i.e. y = In a, it is 
practical to recall, for future convenience, the following relation 



— — 6\iM 

dy 



(/',[,.(! - <). Xi ) + ( - -w de )$>,< 



(2.10) 



which can be used for changing the integration (or derivation) variables from y to Qm itself. 
In the ACDM model the illustration of the analytical and numerical results will be given in 
terms of the following fiducial set of parameters determined on the basis of the WMAP 7yr 
data alone HEJE]: 

(O b0 , ^co^deo, h , n B , e re ) = (0.0449, 0.222, 0.734, 0.710, 0.963, 0.088), (2.11) 

where, by definition, f^Mo = ^co + ^bo; furthermore e re denotes the optical depth at reion- 
ization and n s is the spectral index of curvature perturbations assigned as^ 



97T 2 -» / k \ n s — 1 

(K(k,T)K(p,r)) = —V^d^ik + p), Vn {k) = An(^-J ■ (2.12) 



For the set of parameters of Eq. (12.111) the amplitude of the curvature perturbations at 
the pivot scale k p is given by An = (2.43 ± 0.11) x 10~ 9 . It must be stressed that the 
conventions adopted here to assign the power spectra of curvature perturbations coincide 



7 The specific definition of 1Z in terms of the variables employed in the present paper is reminded hereunder 
at Eqs. (f272l]) and (f^22|) . 
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with the conventions of the WMAP collaboration and will be consistently followed hereunder 
when assigning the power spectra of the fully inhomogeneous magnetic fields (see Eq. (13. 8p 
and discussion therein). If Wd c 7^ — 1 the fiducial set of parameters determined from the 
same observational data (but assuming u>CDM scenario) is given by: 

(O b0 , fi c0 ,^deo, h Q , n s , e re ) = (0.052, 0.26, 0.67, 0.69, 0.952, 0.084), (2.13) 

with Wde = — l-20lo|g an d -An = (2.50 ± 0.14) x 10~ 9 . By analyzing different data sets in 
the light of the tuCDM scenario, the error bars on Wd e either increase or they are restricted 
to an interval pinning down the ACDM value Wde = — 1- In the numerical applications 
and for illustrative purposes we shall often fix the cosmological parameters as in Eq. (12.111) 
by allowing for a variation of Wd e - Values Wd e < — 1 will be excluded for practical reasons 
since, in these cases, the background may evolve towards a singularity in the future [57] 158]. 
This is, of course, only a practical choice since, in principle the present discussion can 
also be extended to the situation where Wde < ~ 1 where interesting models may exist 

[33 EH m\ EH su. 

There are three redshift scales which will be used both in the analytical considerations 
as well as in the numerical discussion: the redshift of last scattering (coventionally denoted 
by the redshift of photon decoupling (conventionally denoted by Zd cc ) and the "drag" 
redshift marking the end of photon drag on the baryons (conventionally denoted by Zdrag)- 
Following the notations already introduced in Eq. (II. ip the value of z* can be expressed as 
(see, e.g. [27J and references therein): 

z* = 1048[l + (1.24x lO- 3 )^ - 73 ^^^^], (2.14) 

0.0783 a; b - °- 238 _ 0.560 

91 " [1 + 39.5 w b ° 763 ]' 92 l + 21.1a; b 1 o 8r { } 

The parameters of Eq. ( 12. lip , once inserted into Eqs. ( I2.14p - (l2.15p imply z* = 1090.77 
in excellent agreement with the direct determination obtained in terms of the WMAP 7yr 
data [U |5j [6] giving z* = 1090.79io;g2. The redshift of photon decoupling is rather close to 
z* so that for the accuracy of the estimates presented here we will effectively consider that 
z* ~ Zdec- For instance the WMAP 7yr data analyzed in terms of the ACDM paradigm imply 
^dec — 1088.2 ± 1.2. The approximate redshift at which photon drag ceases to be effective 
on the baryons can be expressed as 



i oqi .,0.251 

= M0 nS9S [l + &rf ], (2.16) 

1 + 0.659 c^ 28L boJ ' K J 

h = 0.313 w M ° 419 [1 + 0.607 c^ 6 74 ], b 2 = 0.238 c^jf. (2.17) 

Again, the parameters of Eq. (12. lip imply Zdrag = 1020.3 which is compatible with the 
WMAP 7yr data implying Zdrag = 1020.3 ± 1.2. If the evolution of the critical fractions of 
matter, radiation and dark energy is parametrized in terms of y = lna/cide, the parameters 
of Eq. (12. lip lead, with obvious notations, to the following three scales = —6.657, 
Vdec = -6.654 and y diag = -6.589, i.e. approximately y* < y dec < 2/dm g - 
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2.2 Magnetized fluctuations of the geometry 

The only source of large-scale inhomogeneity of the ACDM and wCDM paradigms resides 
in the curvature perturbations whose magnitude determines primarily the normalization of 
the matter power spectrum. The introduction of large-scale magnetic fields in the minimal 
ACDM scenario has been previously scrutinized in a number of analyses either in specific 
gauges (27J I2H1 129] or even within earlier gauge- invariant approaches [43J. The parame- 
ters characterizing the magnetized w;CDM scenario and the ones of the magnetized ACDM 
scenario have been discussed, respectively, in [12J and in [27J. 

The synchronous gauge description is often preferable for the analysis of magnetized 
perturbations as pointed out long ago [29J . Moreover, in the problems related to the growth 
of structures, the synchronous gauge leads to an evolution equation for the total matter 
density contrast which is comparatively simpler than in other gauges. Consistently with 
Refs. [27J |28j [29] the scalar fluctuations of the geometry in the synchronous gauge will be 
expressed, in Fourier space, as: 



5P gij (k,r) = a 2 (r) h(k,r) + 6£(fc,r) - ^) 



(2.18) 



with ki = ki/\k\. Conversely, in the conformally Newtonian gauge the perturbed entries of 
the metric are given by: 

5i cn) goo(k,r) = 2a 2 (T)<P(k,T), 5i cn) g tj (k,r) = 2 a 2 (r)^(A;, r)S tj . (2.19) 

The parametrizations of Eqs. (I2.18P and (I2.19P are related by the appropriate coordinate 
transformations, i.e. 

^(k,r) = -Z(k,T) + ^[h(k,T)+6Z(k,T)]', 

0(fc,r) = -^[h(k,T) + 6ak,T)}" + H[h(k,T)+6ak,r)}'y (2.20) 

The curvature perturbations on comoving orthogonal hypersurfaces, i.e. TZ(k, r), are defined, 
in the gauge of Eq. ( 12. 18ft . as: 

*<*• r) = « + ^ - *<-><*, r) = - (2.21) 

where the arrow denotes the resulting expression of TZ(k, r) in the conformally Newtonian 
gauge obtainable by shifting metric fluctuations from one coordinate system to the other; 
the gauge-invariance of 1Z(k, r) implies, obviously, that lZ(k, r) = TZ^ cn \k, r). From the last 
equality in Eq. (I2.2ip it is also apparent that when TV = 0, = TZ(k). In more general 
terms, by solving Eq. ( 12.211) in terms of £, it can be shown, after integration by part, that: 

ak, r) = TZ(k, r)-^r ^l^'C*. n) dr u (2.22) 
a{T) Jo ri[Ti) 
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where T\ is an integration variable and where the prime denotes, as usual, a derivation with 
respect to r. Both in analytical and numerical calculations the normalization of the curvature 
perturbations is customarily expressed in terms of 71; for this reason Eqs. (I2.2ip and (12.221) 
turn out to be particularly useful in the explicit estimates. The same transformations of 
Eq. (12.201) can be used to gauge-transform the governing equations from the conformally 
Newtonian frame to the synchronous coordinate systeno As a general comment on the 
conventions employed in Eqs. (I2.18P and (12.191) we recall that they are consistent with the 
ones of [1~2||2~T] and they differ slightly from the ones of Ref. [H] (see also [15]). Furthermore, 
in the present approach the curvature perturbations on comoving orthogonal hypersurfaces 
(i.e. 1Z) are often related, in the large-scale limit, to the curvature perturbations on uniform 
density hypersurfaces (see, e.g. [IB] and also [21 E])- Finally, general analyses of large-scale 
magnetic fields performed within the covariant formalism are available in the literature (see, 
e.g. [S]) but they are not central to the present discussion. 

The Hamiltonian and momentum constraints stemming, respectively, from the (00) and 
(Oi) (perturbed) Einstein equations are given, in real space, by 



2V 2 £ + Hti = SirGa 2 



V 2 £' = ^Ga l 



(pt + pt)0t + 



4pf + 5 s p de + 5 s p B + 5 s p E 
Ex B 



4vra 4 



(2.23) 
(2.24) 



where, 5 s pf = 5 s p^ + 5 s p c + 5 s p u + 5 s p 7 denotes the density fluctuation of the fluid sources 
in the synchronous gauge. In Eq. ( I2.24p the three-divergence of the total velocity field, i.e. 
9t = V • Vt can be expressed in terms of the contribution of each individual specie^] 

4 4 

(Pt + Pt)#t = gPA + gP 7 # 7 + PiA + PA + (w dc + l)p dc 9 de . (2.25) 

In Eqs. (12.231) and (12.241) the effects of the large-scale electromagnetic fields have been 
included in terms of the comoving electric and magnetic fields E(x,r) = ci 2 (t)£(x,t) and 
B(x,t) = a 2 (T)B(x,T), i.e. 



B 2 E 2 5 s p B . S s p E 

S7ra 4 87ra 4 3 3 



where B 2 = \B\ 2 and E 2 = \E\ 2 . For z > z dvag the baryon and photon velocities effectively 
coincide with 9^ 

b+^k , b = !^i±M , (2.27) 

m e + rrii 



8 Not only the metric fluctuations will change under coordinate transformations but also the inhomo- 
geneities of the sources. In particular, it can be easily shown that S^p = S^p - p'[(h' + 6£')/(2fc 2 )]. These 
considerations are relevant in connection with the last part of section [3J 

9 Following exactly the same conventions established in Eqs. (|2.2|) and (|2.3j) . the various subscripts denote 
the velocities of the different fluid components; recall that, for a generic species X, the notation 9x — V -vx- 
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and Eq. (I2.25P can be written as: 

(Pt + fh)Ot = tpuO V + + Rb)6yb + PA + (™de + l)Pde0de, #b(«) 



3 3 4w 7 o(Zde + 1) 

(2.28) 

i?b(«) denotes the well known weighted ratio between the baryonic matter density and the 
photon energy density determining, for instance, the acoustic oscillations in the temperature 
autocorrelations of the CMB. The evolution equations for E and B are: 

V • B = 0, V • E = 4vrp q , (2.29) 
VxE + B' = 0, VxB = 4ttJ + E' (2.30) 
p q + V-J = 0, p q = e(rii - n e ) , J = e(riiVi - n c v c ). ( 2 -31) 

The total current J has been expressed in terms of the two-fluid variables, however, as in 
the case of flat-space magnetohydrodynamics, J can be related to the electromagnetic fields 
by means of the generalized Ohm law obtained by subtracting the evolution equations for 
the ions and for the electrons 



J = *(E + tbXB + ^-—), a=— , „ (2.32) 

V en n e J 47r{ar ie + (4/3)[p 7 /(n m e )jr e7 } 

where a is the conductivity; r ie and T 7e are, respectively, the electron-ion and electron- 
photon interaction rates. The three terms appearing in the Ohm's law are, besides the 
electric field, the drift term (i.e. x B), the thermoelectric term (containing the gradient of 
the electron pressurd^l) and the Hall term (i.e. JxB). For frequencies much smaller than 
the (electron) plasma frequency and for typical length-scales much larger than the Debye 
screening length the Hall and thermoelectric terms are subleading for the purposes of the 
present analysis both before and after photon decoupling [27] . 

In the resistive magnetohydrodynamical description adopted in the present paper the 
electric fields are subleading in comparison with the magnetic fields by a dimensionless factor 
going approximately (L o") -1 where L is the typical length-scale which can be identified with 
the wavelength of the fluctuation. For z < z^ec the ionization fraction drops and by Zdrag the 
conductivity can be estimated as 

^4.36xlO-evf^V / Y^V /2 f^±iy /2 , (2.33) 
V2.725K7 V0. 13347 \z cq + l) ' v ; 

which is a pretty large value as it can be argued by computing explicitly the suppression 
parameter arising, for instance, in the magnetic diffusivity terms 

P / h \ 2 / 7i \ -V 2 / hr 

o 1/^—20/ \ I drag \ / "C 



3xHT 2U -— -^^) — , (2-34) 
Tia \MpcJ V1020.3/ V0.71,' 



ln In Eq. (|2.32|) the comoving electron pressure is given by p c = n c T c where n c and T e are, respectively, 
the comoving concentration and the comoving temperature of the electrons. 
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where we estimated ^M^drag) — 0.350f2Mo- All the terms explicitly suppressed by the 
conductivity and by one or two spatial gradients of the magnetic field will therefore be 
neglected since they are subleading in comparison with the other terms. In spite of the 
smallness of the ionization fraction after decoupling (and after the drag time) the plasma 
approximation is still very good and the plasma is an excellent conductor as Eqs. (I2.33p - 
(I2.34p clearly show. When the ionization fraction x e drops almost suddenly from 1 to about 
lO^WMo/^bo [MIHS] (see also [501 EE] f° r a delayed recombination scenario) the concentration 
of the free charge carriers diminishes. Since the free charge carriers drop faster than the 
temperature the Debye scale increases. Overall, however, the plasma parameter decreases 
since 



plasma = ~T? ~ = ^ J^-y/x^ = 2.3 X 10"^f ^ (2.35) 

VDn x e V 7T \0. 02258/ 

where Vd = 47rA D /3 is the volume of the Debye sphere, Ad is the Debye length and £(3) = 
1.202. Equation (I2.35P shows that the plasma approximation is still rather accurate for 
z < -^drag! furthermore the largeness of the conductivity justifies the approximations adopted 
in the present analysis. At least in the general equations written below the terms involving 
the electric fields will be kept both for sake of completeness and for future convenience. 

The electromagnetic pressure and all the sources of anisotropic stress enter the perturbed 
(ij) components of the Einstein equations: 

h" + 2Uti + 2V 2 £ = 2AnGa 2 [5 sPi + 5 sPde + 5 s p B + 5p E } , (2.36) 
1 



^[a 2 (h + 60' ]' + 2V 2 £ = 327rGa 2 



a 



3 / 

PvO v + p 7 CT B + P 7 OE + ^Pdel^de + l)0"de 



,(2.37) 



where 5 s pf = (5 s p 7 + S s p v )/3, in analogy with S s pf, denotes the fluctuation of the pressure of 
the fluid components. In Eq. (I2.37|) the notation 

didjUf = ^p u V 2 a u + ^p 7 W B + ^p 7 V 2 a E + p de (^de + l)V 2 a de , (2.38) 

has been adopted where, as usual, the total anisotropic stress 11^ has been made explicit in 
terms of the anisotropic stresses of the different species: 

u{ t = ul + n^ B + n| E + u{ de . (2.39) 

The various subscripts in Eq. (12.391) denote the corresponding components and, in particular 
the electromagnetic contribution: 



W — 



j 



BiB 3 - -^B 2 
3 



IT -J- 



j 



8 3 - 

EiE 3 - —E 2 
3 



(2.40) 



The species of the plasma either interact strongly with the plasma (like the elctrons, the ions 
and the photons) or they only feel the effects of the geometry (like the CDM component and 
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the dark-energy). The baryon evolution equations are given by: 



^ = % - K o', + ne' h = V ' [ ^ B] + \^{e, - <h>), (2. li : 

l a p b -J Pb 



while the evolution equations for the photons are: 

= -~V 2 5 7 + e'(6 h - ft,), 5; = - ^6 T (2.42) 

In Eq. (I2.42p e' = anocr e7 is the differential optical depth due to electron-photon scattering. 
The CDM density contrast and velocity will evolve, respectively, as 

5' C = ^-9 C , 6' C + H6 C = 0, (2.43) 
while the evolution of the neutrinos is given by: 

^ = -Iv 2 <f„ + W 5' v = lh'-^0 v , (2.44) 

4 3 2 4 

<r' = —0 V - — b! - (2.45) 

y 15 10 3 15 5 s ' v ; 

where JF V ^ reminds of the coupling of the monopole and of the dipole to the higher multipoles 
of the neutrino phase space distribution. The latter term will be set to zero in the class of 
initial conditions discussed in the present paper but it can be relevant when magnetized 
non-adiabatic modes are consistently included. 

As far as the dark-energy fluctuations are concerned, the situation changes qualitatively 
between the ACDM scenario and the u>CDM case. If w dc = — 1, 5 dc = ftic — and the dark 
energy does not cluster on subhorizon scales. If w dc ^ — 1 the dark energy perturbations 
evolve; to prevent instabilities the barotropic index w de and the sound speed c de are assigned 
indipendently. The latter choice implies that the total pressure fluctuation inherits a non- 
adiabatic contribution which is proportional to (w de — c^ e ); thus the fluctuations of the dark 
energy pressure can be written as [52] (see also {53"| 154]): 

Q 

S s Pde = cl e 6p de + 5jW, 5p nad = 3H(1 + W dc ) (c\ c - W dc )p de -j, (2.46) 

from which it is clear that c\ e is the sound speed in the frame comoving with the dark energy 
fluid. Using Eq. (12.461) together with the (perturbed) covariant conservation equation for the 
energy-momentum tensor of the dark energy the following pair of equations can be readily 
obtained in Fourier space: 

+ 3^(4 - w de )S dc + (w de + 1) j [fc 2 + m 2 (c 2 de - w de )\ ^ - ^} = 0, (2.47) 
d' dc + Ui\ - 34)fta e - = 0, (2.48) 
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In Eq. (I2.47P the possible presence of an anisotropic stress for the dark energy contribution 
has been neglected. The bounds on c de are currently rather loose and we shall assume, as 
customarily done, that < c de < 1. The expression reported in Eq. (I2.46[) is mathematically 
analog to (but physically different from) the scalar fluctuation of the total fluid pressure 
which can be written as 

8 s Pi = c 2 st 5 sPi + 5p nadJ c 2 st = % (2.49) 

Pt 

where c 2 t is now the total sound speed of the ordinary fluid sources (i.e. characterized by 
a barotropic index which is positive semidefinite) and 5p nad parametrizes the non-adiabatic 
pressure fluctuations stemming from the spatial variation of the chemical composition of the 
plasma. In the standard case 5p nad ^ for the CDM-radiation mode, for the neutrino-density 
mode, for the neutrino velocity mode and for the baryon-radiation mode. All these modes 
can be appropriately generalized to include the contribution of fully inhomogeneous magnetic 
fields [3B], however, for the present purposes we shall avoid the technical complication of the 
non-adiabatic modes and assume throughout adiabatic initial conditions (i.e. £p nad = 0) both 
at the level of the Einstein-Boltzmann hierarchy and at the level of the growth equation. 



3 Estimates of the magnetized growth rate 

The evolution equation for the growth rate introduced in Eq. (11. ip will now be derived and 
solved in different approximations. The combination of Eqs. (I2.4ip and (I2.43P leads to the 
following equation for the matter density contrast: 



5" + U5' 



wmo a 4 Pb 2 



(3.1) 



In (13.11) no approximations have been made on the relative weight of the relativistic correc- 
tions. By summing up Eqs. f !2.23|) and (I2.36P the following equation can be obtained 



h" + Uti = 8nGc 



5 s pi + 35 s p f + 5 s p dc + 35 s p de + 5 s p B + 35 s p B 



(3.2) 



The pressure fluctuations of the fluid sources have been separated, according to Eq. (12. 49 p . 
into the adiabatic and non-adiabatic contributions; the non-adiabatic contribution of the 
fluid sources will be set to zero and, therefore, inserting Eq. (13 .2p into Eq. (13. ip . the result 
is: 



5'' + V.5' 



WboV-(JxS) 3 2 r Q R Q r 

+ - ft \ly{ \ m + 2K^-—\Ib + 7^—0 



+ 



a 4 p b 



M 



R 



M 



<5de(l + 34) + 9U(1 + w dc )(c 2 de - w de )^ 



(3.3) 
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At the right hand side of Eq. (I3.3P the magnetic energy density is expressed in units of the 
photon energy density by defining S s Pb = ^b(x, r)p 7 . Equation (13. 3 j) can be further modified 
by recalling that, on the basis of simple vector identities (see first paper in [56J) 

3 V-(Jx£) _ didjtfj v 2 n B 
4 a 4 p 7 P 7 + P 7 4 

For a 3> a cq , f2 R <C fi M and f2 R ^de in Eq. (13. 3p . Moreover, assuming |<5 dc | <C \S m \ the 
approximate form of Eq. (13. 3ft becomes: 

C + - t^iA. = - Uh \ [4VV B - V 2 fi B ]- (3.5) 

Equation (13. 5ft only applies for redshifts z < Zdrag and will now be used for analytic estimates. 
In the limit Qb — > and <tb — > when Eq. (13. 5p reproduces the standard equation analyzed 
in various situations for the calculation of the growth index in diverse models of dark energy 
[1"3~| US] (see also [581 EH EHl EJJ |62]). By changing the variable from the conformal time 
coordinate to the natural logarithm of the normalized scale factor, Eq. (13. 5 p becomes, for 

Vi ^* Z/drag) 



d 2 S m 
dy 2 



IT - l Wm = 'J ^TWn^B - V 2 ^b], (3.6) 



where, recalling the expressions of Eqs. (I2.5P and (I2.28p . 



r> / \ 3 WbO 



4 w 70 z dc + 1 

2^de + 1 



W(y) = Ut^ L —n uo e 



-y 



[z eq + 1) + (^de + ^e-^ + (z ccl + l)e- 3w ^ 



(3.7) 



Since the magnetic fields are stochastically distributed, the ensemble average of their Fourier 
modes obeys: 

-» 2lT 2 -> f k \ n B — i 

(BiWBjtf)) = —P^VBi^ik + p), V B (k) = AB[j-J , (3.8) 

where Pij(k) = (k 2 5ij — kikj)/k 2 ; Ab the spectral amplitude of the magnetic field at the 
pivot scale k^ = IMpc -1 . The conventions adopted in Eq. (13. 8 p are the same as the ones 
employed in Eq. (I2.12p reproduces the conventions followed by the WMAP collaboration 
(see, e.g. [H El [6] and earlier data releases). There are some who assign the power spectra of 
curvature perturbations as in Eq. (" 12 . 121) (where the scale-invariant limit is n s — > 1) while, on 
the contrary, the power spectra of magnetic fields are assigned in such a way that their scale- 
invariant limit would correspond to tib — > —3 (and not to ub — >■ 1). It seems deliberately 
confusing to use different definitions for the same mathematical object in the same physical 
framework. In other words, denoting with n a generic spectral index, it is certainly possible 
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to assign the power spectra appearing in Eqs. (I2.12p and (13.81) by engineering the scale- 
invariant limit for n — > —3 rather than for n — > 1. However this choice must be consistently 
implemented for all the power spectra. In the present paper the power spectra of Eqs. 
(I2.12p and (13. 8p are then assigned with the same conventions. Consequently, in Eq. (13. 8p . 
Ab has the correct dimensions of an energy density!"! Furthermore Ab can be related to the 
regularized magnetic field intensity B^ [T2"l |2"7] which is customarily employed to phrase the 
comoving values of the magnetic filed intensity. In the case when tib > 1 (i.e. blue magnetic 
field spectra), Ab = (2n) nB ^ 1 Bl/Y[{riB — l)/2]; if n B < 1 (i-e. red magnetic field spectra), 
Ab = [(1 — ^b) /2](k p / k^)^ nB ^ Bl. In the case of white spectra (i.e. Ub = 1) the two-point 
function is logarithmically divergent in real space and this is fully analog to what happens 
in Eq. (I2.12p when n s = 1, i.e. the Harrison-Zeldovich (scale-invariant) spectrum. 

The parameter space of the magnetized wCDM models and of the magnetized ACDM 
scenario have been investigated, respectively, in [12] and [27] . In a frequentist approach, 
the boundaries of the confidence regions obtained in [TJJ [27] represent exclusion plots at 
68.3% and 95.4% confidence level. When moving from the magnetized ACDM scenario to 
the magnetized wCDM model we have that the the parameters maximizing the likelihood 
get shifted to slightly larger values^ 

(tib, £l)acdm = (1-598, 3.156nG) ->■ (n B , B l ) wC bm = (1-883, 4.982 nG), (3.9) 
(n B , 5 L ) ACDM = (1.616, 3.218nG) -> (n B , B L ) wCDM = (1.913, 5.163 nG). (3.10) 

Even if the addition of a fluctuating dark energy background pins down systematically larger 
values of the magnetic field parameters, the results of [121 EZ] will be used here just for a 
consistent illustration of the results. If n B > 1 (as in Eqs. (I3.9p - (l3.10p ) a useful analytical 
approximation of the power spectra of <jB{k,y) and QB{k,y) can be written as 

(n B (k,y)Mp,y)) = ^-MKy)^\k+p), 

2tt 2 

(aB(k,y)a B (p,y)) = —V a (k,y)6 (3 \k + p), (3.11) 
where Qbl = Bl/(8np^) and so the power spectra are given by: 

Vn(k,y) = T(nB,k D )n BL ^—j , V a (k, y) = G(n B , k D )n BL ^— j . (3.12) 

The functions J-"(rzB, ^d) and Q(n B , k D ) are defined in terms of the spectral index ub and of 
the diffusive wavenumber k^: 

4(7-n B ) 8 /fc\^ n B 4 / k \ 5 ~ 2 "b 

(nB ' D) ~ 3{n B - 1)(5 - 2n B ) ~ 3{n B - 1) \k~J + 2n B - 5 \k^J 



11 In other words, as observed long ago in a closely related context [55] . the magnetic power spectra assigned 
here coincide with the magnetic energy density per logarithmic interval of frequency. 

12 The difference between Eqs. (|3.9|) and p.lO[) is that the parameters of Eq. Q3.9p are obtained from the 
analysis of the temperature autocorrelations while the parameters of Eq. (|3.10|) are obtained by adding the 
data points of the cross-correlations between temperature and E-mode polarization [12 . 
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G(n B ,k D ) 



n B + 29 



1— ris 



7 



5(2n B -5)U D 



5— 2?ib 

(3.13) 



15(5 — 2n B )(^B — 1) 3(n B — ' 
Recalling that the angular diameter distance at last scattering can be written as D^z*) = 
2d\(z*)/ (H \/Q M0 ), k A (z*) and k-£>(z*) are determined in terms of the parameters of the 
fiducial set of parameters (given, for instance, by Eqs. (12. lip and (I2.13P ): 

MaO 2240 d A CO /z,\5/4 



k A (z*) = l/D A (z*), 



'?"r* 



103 



0.24, ,-0.11 
M 



UJ W UJ 



(3.14) 



The results of Eqs. ( 13. 12ft and ( 13.131) compare pretty well with the numerical solution of 
the two convolutions defining Vn(k,y) and V a (k,y) [29J, |56] . Note that, incidentally, the 
same kind of integrals determining Vn{k,y) and V a (k,y) arise when computing secondary 
graviton spectra from waterfall fields where the recent numerical results of [57] show, once 
more, excellent agreement with the semi-analytical scheme leading to the results ( 13.121) and 
(I3.13p . The relative magnitude of Qbl and of An (i.e. the amplitude of the power spectrum of 
curvature perturbations) depends on B^ (i.e. the typical magnetic field intensity regularized 
over a typical length-scale k^ 1 ) 



n 



BL 



A 



39.56 



n 



nG 



T- 



70 



2.725 K 



A 



n 



2.41 x 10" 9 



(3.15) 



In Fourier space and with all the specifications given above Eq. f]3. 61) can be written as: 



d 2 5 n 



where 



dy 2 
S(k,y) 



1 3 n 

2 - 2 Wdc(1 



n 



M / 



d5 m 
dy 



k 2 u h0 



S{k,y), 

4 



(3.16) 



(3.17) 



R h (y)H 2 (y)uj M o - 

By introducing the explicit definition of the growth rate Eq. ( I3.16P and ( I3.17P lead to an 
integro-differential equation whose explicit form can be written as: 



dy 



- - 7;Wde(l - ^ 



/ - 



1 



■VL M U[f] y h y), 



(3.18) 



2W B /3 

where U[f] y- u y] is a functional of the growth rate; W B = W B (k, n B , Bi) depends upon the 
magnetic field spectra and upon the spectrum of matter inhomogeneities at y\\ 



U[f\ y h y] = t 
W B (k,n B ,B L ) 



-I(Vi,V) 



i(yu y) 



f(x) dx, 



3w mo \5m(k,yi)\ \ Hi 



k he 



2 r 



A\a B {k)\-\n B {k)\ 



(3.19) 
(3.20) 



For typical scales k 3> k eq the matter power spectrum at y eq can be approximated a 



P s (k,y eq ) = T 2 (k/k eq )5 2 n (k) -> -An 



n s — 1 



ln 2 (fc/fc eq ), 



(3.21) 



13 In Eq. (|3.2 1[) Jfl-(fc) denotes the initial spectrum which is related to the spectrum of curvature pertur- 
bations while T(fc/fc cq ) denotes the transfer function. 
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so that Eq. (I3.20p can also be written as: 

5u 7 o / k \ nB_1 f k \ i 1 -™*)/ 2 n BL 



W B (k,n B ,B L ) 
jO B {n B , k-o) 



6 u'mq V k L 



H$ In (fc/fc, 



£b(^b, fco) 



(2tt 



I TIB— 1 



r[(n B - l)/2] 



4Vl^(n B ,A; D ) 



(3.22) 



The evolution variable in Eq. (13. 18ft is y but since the relation between Qu and y is regular 
and invertible, Eq. ( 12 .101) can be used to obtain the evolution of / directly in terms of Qu' 



3u>hp(1 — f2 



df 



M, 



where 



I (to 



c/ln^M 
1 



+ / 2 + 



1 3 

2 ~ 2 

/(«) 



W de (l - f) 



/--fi M = W B W[/;fi M ],(3.23) 



^ B = ^ B (fc,n B)J B L ) 



3w d e ^Mfei) ^(l - 0) 

The solution of Eq. (I3.18P and f !3.23|) will now be parametrized as 



/(*,!/) 



f(y) 



l-2W B U[f;y h y}/3' 



(3.24) 



(3.25) 



where f(y) obeys the following equation: 



dy 



- - 7^ dc (l - Q M ) 



f ~ = 0, 



(3.26) 



and where U[f ; yi, y] has been defined in Eq. (13.191) : notice, however, that in Eq. (13. 25ft U 
is a functional of / and not simply of /. To prove that Eqs. ( 13.25j) and (I3.26P indeed solve 
Eq. (I3.18P and f!3.23j) let us start by showing that if Eq. (I3.25P holds, then the following 
relation is also verified: 



U[f\y h y] 



(l-2W B /3)U[f ]yu y] 
\-2W B U[f-y h y}/3 



(3.27) 



Equation (13.271) follows by integrating Eq. (13.251) over a dummy variable (be it x) between 
y\ and y; let us then show this explicitly and first rewrite Eq. ( 13.25ft as 



f(k,x) = f{x) + 



2W B V[f;x il x]f(x)/3 , 
l-2W B U[f;x h x]/3 : 



(3.28) 



note that, at the right hand side of Eq. (13.281) . the argument of U is / (and not /). By 
integrating over x the left and right hand sides of Eq. (I3.28j) between y^ and y the following 
equation can be easily obtained: 



f(k, x) dx 



f(x) dx — In 
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l-2W B U[f-y h y}/3 
1 - 2 W B /3 



(3.29) 



where the logarithmic contribution arises by direct integration of the second term in Eq. 
(I3.28P between the two limits y- x and y. The left and the right hand sides of Eq. (I3.29j) can 
then be exponentiated and Eq. (13.271) is recovered. Note that, according to Eq. (I3.27p . 
U[f ; yi] = 1 as it must be for consistency with the definition of U[f ; yi, y] given in Eq. 
(I3.19p . Bearing now in mind the results of Eq. (I3.27p . f(k,y) given by Eq. (13.251) can be 
plugged into Eq. (13.181) and the obtained equation for f(y) turns out to be exactly the one 
given in Eq. (13.261) . 

To derive an explicit solution for the growth rate we need to solve Eq. (I3.26p . but this 
part of the problem is more conventional. By assuming that / = Q]^, Eq. (13.261) becomes: 



3w de (l — ^ 



In ft 



c?7 



dlnft 



+ 7 



M 



o 7 



2 " 2 Wdc 



1 - ft 



o 7 



3 ft 



0. (3.30) 



Eq. (I3.30p can be solved as a power series in e = 1 — Qm by neglecting, to lowest order in 
e, the derivatives of 7 with respect to ftM- A posteriori this assumption will turn out to be 
rather accurate and it is commonly adopted in the absence of magnetic field contribution 
[UJdS]. Thus, from Eq. (13.301) the following relation can be determined: 



3w de (l - ft M )7 + + 



-w de (l - ft 



ft 1 " 7 



0. 



By expanding the left hand side we obtain, to leading order in e = 1 — ft 



7 



3w de -3 3 (1 - w de )(l - 3w de /2) 



Qw 



+ 



dc 



125 



(1 - 6w de /5) ; 



£>(e 2 



(3.31) 



(3.32) 



Inserting the leading order result of Eq. ( 13.321) into Eq. ( 13.25}) the explicit form of the growth 
rate becomes 



f(k,y) 



o 7 

M 



2W B U[f;y u y]/3' 



u[f;yu y] = exp 



^il A (x)dx 



(3.33) 



where, now, 7 = (3w de — 3)/ (6w de — 5) + 0(e) and, as already mentioned prior to Eq. (I3.20p . 
Wb = W / b(^, ""-b, -Bl) is given by Eq. ( 13. 22H . Even if more accurate analytical results can 
be obtained by including higher orders in e, to privilege the simplicity of the approach and 
of the comparison we will stick, for the present treatment, to the explicit expression of Eq. 
(I3.33P with the growth index evaluated to lowest order in e. The integral appearing in Eq. 
(I3.33P can be performed either directly or by using Q as integration variable (see, e.g. Eq. 
(I2.10p for the appropriate change of integration variables): 



(x) dx 



ft 7 



3w de JQ M (yi) ft (ft - I] 



dQ. 



(3.34) 



The lower limit of integration indicated in Eq. ( I3.30|) can be estimated by recalling that 

I- ^RO 



cqj 



2H 


ftRO 


f ftdeO ^ 3w de 


-1 




ftMO 


VftMO/ 





^2' 



ft 



.MO 



c(io- 5 ). 



(3.35) 
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Thus, if y eq C ?/i < ya e , ^u(yi) > 1/2- For instance y cq = —7.78 for the fiducial set of 
parameters of Eq. (12. lip . Thus y\ can be safely chosen in the range —7 < y\ < —4 if we 
want to be consistent with the approximations made so far. In Fig. [1] the result of the 
numerical evaluation of Eq. (I3.34p is reported in the two physical cases which have been 
taken as extreme, i.e. w<ie = —1 and w^e = —0.3. With the continuous line the analytical 
interpolation /(?/;, y) = (y— ?/drag) is reported. The fiducial set of the cosmological parameters 




Figure 1: The result of the numerical evaluation of the integral of Eq. (13.341) is compared 
with the analytical interpolation discussed in the text (plot at the left). The agreement is 
reasonable for most of the ranges of Wde and y but it gets worse (see the right plot) in the 
region of large y (i.e. y — > 1) and large w& e (i.e. w& e — > —0.3). 

coincides with the one reported in Eq. (12. lip . The values of and n-Q lie within the range 
of parameters determined from the analysis temperature and polarization anisotropies of 
the CMB within the magnetized ACDM scenario and in the light of the WMAP data alone 
[13 [271 ( see also Ec i s - (BSD and ( ETTDjl ). 

In Fig. [2] (plot at the left) the growth factor of Eqs. fl3T25|) and (13331) - f[3T34p is illustrated. 
The thick lines in Fig. [5] correspond to the numerical evaluation of the integral of Eq. (I3.34p 
while the thin lines are obtained by means of the analytic approximation (full line in Fig. 
[p. The analytic approximation is adequate for the purposes of the present analysis. In the 
right plot in Fig. [2] the common logarithm of the growth factor is illustrated for a fixed scale 
and fixed magnetic field amplitude but different spectral indices. In the left plot of Fig. [3] 
the growth rate Eqs. (13.251) and (I3.33I) - (I3.34I) is plotted for different values of Wd e and by 
fixing all the other parameters to a selected fiducial value mentioned in the legends and in the 
labels appearing at the top of each plot. The range of y has been narrowed in the plot so that 
the curves corresponding to different values of w^e are more clear. In the right plot of Fig. 
|3]the direct numerical solution of Eq. (13 .6p is illustrated with a full (thin) line. The dashed 
and dot-dashed lines denote, respectively, the semi-analytic and the fully analytic result. 
By semi-analytic result we mean the expression for the growth rate obtained by estimating 
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Figure 2: The growth factor is illustrated analytically and numerically (plot at the left) for 
different values of the comoving wavenumber and of the magnetic field intensity but keeping 
the magnetic spectral index ub fixed. In the right plot the common logarithm of the growth 
rate is illustrated for different values of the magnetic spectral indices by fixing all the other 
parameters, as explained in the text. 



A = 2.43x 1 0" 9 , n = 0.963, B. = 1, n„ = 1.3, k =0.3 Mpc" 1 A = 2.43x 1 0" 9 , n =0.963, w . = -1 , B, = 1 nG, n=1 .5, k=0.03 Mpc" 1 

H S L b H S Ge L d 




-0.4 -0.2 0.2 0.4 0.6 0.8 1 -5 -4 -3 -2 -1 1 



y y 

Figure 3: The growth factor for different values of w^ e (plot at the left); the full numerical 
solution of the evolution equation of the density contrast is compared with the analytic and 
semi-analytic approximations previously discussed (plot at the right). 



numerically the integral of Eqs. (I3.25P and (I3.34p (see also dot-dashed and dashed lines in 
FigJT]). in the analytic case the integral I(y, y\) is estimated by the approximation illustrated, 
with the full line, in FigJTJ The direct numerical solution of Eq. (I3.6P leads to a result which 
is approximately located between the analytical and semi-analytical approximations. Other 
choices of the relevant parameters lead to results whose accuracy is comparable with the 
encouraging results of FigfSJ The only caveat with the latter statement is that, of course, 
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Eqs. (I3.25P and f)3.33p -f l3.34p are in good agreement with the numerical solution of the 
approximate equation (i.e. Eq. (13. 6p ). This aspect will be further deepened in section [U 

The results reported so far have been derived within the synchronous gauge description. 
So, for instance, in a different gauge Eq. (13. 5p will assume a different form since 5 m is not 
gauge invariant. This is not a problem since once the calculation is performed in a given 
gauge, the results can be translated in a different coordinate system. Still it is relevant to 
point out that the growth equation obtained in Eq. (13.51) is not gauge-invariant since the 
matter density contrast 5 m does change for infinitesimal coordinate transformations. Very 
often, indeed, Eq. ( 13. 5p is derived in the conformally Newtonian gauge but this demands, 
as we shall show, different assumptions on the relative smallness of the relativistic correc- 
tions. So, even if this is a common problem also in the conventional case the considerations 
developed hereunder seem to be appropriate. The same steps leading to Eq. (13.51) can be 
repeated in the longitudinal gauge with the resuhQ 

K + UC = 3« + W) + V 2 - ^- ( f Xjg) , (3.36) 

^MO a Pb 

where 5 m now denotes the density contrast in the longitudinal gauge; the definition of (f) and 
ip can be found in Eq. (I2.19p . Using the Hamiltonian constraint in the longitudinal gauge, 

^M^m + ^R^R + ^dc^de + -R 7 ^R^B > 



V> = + K0) + -H 

as well as the dynamical equation for ip, Eq. ( 13.361) becomes 

~C + Ht m = -3(W + 2U'ip) - ^-(Jx£) 

wmo a Pb 

3 



(3.37) 



g Q, _ 

2Q R 5 R + fi M 5m + 2R^Q R Q B + Q de (l + 3c 2 dc ) - ^-yWde(Wde + 1)0, 



do 



,(3.38) 



which coincides with Eq. ( 13. 5 j) only in the non-relativistic limit (i.e. (V 2 '?/'! 3> ip" and 
(V 2 ^! 3> *Hip') and in the case ip — > <fi. It could be naively expected that the same physical 
approximations leading to the growth equation in one gauge would lead to the same growth 
equation in another gauge. This naive expectation is incorrect as Eqs. ( 13 .5p and ( I3.38P 
show. The correct conclusion drawn from the comparison of Eqs. ( 13 .5p and f!3.38j) is that 
the non-relativistic limit is implemented in different ways in different coordinate systems. 
The equation for <5 m , valid in the synchronous gauge, holds under milder assumptions, i.e. 

<^ Vlyi and 5 dc — 0. The equation for 5 m , under the same assumptions, also contains extra 
terms which can only be neglected by enforcing the non-relativistic limit in Eq. (I3.38P and 
by consequently neglecting "^CHip' + 21-L'ip). The synchronous description seems therefore 
more appropriate for the consistent computation of the growth rate. It might seem that 
the problem of the ambiguity in the definition of the growth equation could be solved by 
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Note that Sx denotes the density contrast of the species X in the longitudinal gauge. 
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appealing to the standard gauge-invariant descriptions. Consider for instance a conventional 
set of gauge-invariant variables such as the standard gauge-invariant generalization of the 
longitudinal gauge variables i.e. ip — >■ ty, <fi — >■ $ and 5x — > Sx where $ and \I/ are the two 
Bardeen potentials. In this case the equation for the gauge-invariant density contrast 6ffl 
will be the same as Eq. ( I3.38p . It seems more interesting to introduce the gauge-invariant 
variables (see e.g. [SB], first paper) 

Cc = -^ + |, Cb = -^ + y, (3.39) 

which are essentially the density contrasts of the the CDM and of the baryons but on the 
hypersurface where the curvature in unperturbed (see, e.g. (2J [3]). Recalling, from the 
Hamiltonian constraint, that V 2 ip = 12irGa 2 (pt + Pt)(C ~ 7£) we have that the evolution 
equation for ( m = (u; c0 /u;mo)Cc + (^bo/^MoKb becomes 

C + nC = (H 2 - U')(( -K)- 3^ V " ( ^ Xi?) , C = E -Ca, (3.40) 

wmo a Pb a Pt> 

where 1Z is the standard variable describing the curvature perturbations on in the comoving 
orthogonal gauge and ( is the total density contrast in the uniform curvature gauge. This 
description can be used for the computation of the growth rate with some advantages which 
are, however, not central to the present discussion. 

Before concluding this section it is useful to is useful to introduce the magnetic Jeans 
length which allows to express the normalizations of the magnetic power spectra. The 
comoving magnetic Jeans length can be defined as 



A BJ = cL^, cl=J§r, (3.41) 

where p b = a 3 pb simply denote the comoving baryonic density. From Eq. ( 13.41 j) the explicit 
form of the square of the magnetic jeans length, becomes 



or, even more explicitly, 



Abj = 1.90x10- 2 (^^) (^)Mpc. (3.43) 
J V0.02258/ VnG/ P V ; 

In terms of the magnetic Jeans wavenumber fc B j = 2vt/Abj the quantity M^(Mb)5l) can 
be written as 

W B (k,n B ,B L ) = f * Vf f *V'- W ^»^°> , (3.44) 



WMoVfcBj/ \k L J \k p J ^/Ar, hi (k/k e(l ) 
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The expression (13.441) generalizes to the ACDM situation the analog expressions arising in 
the absence of a dominant adiabatic mode seeding structure formation. Equation (I3.44p 
shows that the effect of scale variation of the growth factor is partially dominated by scales 
close to the magnetic Jeans scale at least in the case n s ~ 1 and 1 < tib < 1.5, which is what 
the plots of Figs. [2] and [3] show. 



4 Numerical discussion 

The results of section [3] neglect the contribution of the radiation and of the relativistic 
corrections; furthermore they also assume that the evolution of the dark energy fluctuations 
does not appreciably contribute to the final shape of the growth factor. In the absence of 
large-scale magnetic fields, the first two assumptions are rather reasonable for low redshifts 
and for modes which crossed inside the Hubble radius prior to equality, the third assumption 
is quantitatively correct only in ACDM context but not in the wCDM case. If the plasma is 
magnetized the system changes qualitatively and quantitatively both before and after photon 
decoupling. The numerical analysis of the present section is intended to complement and 
corroborate the results of the previous sections. 

The numerical integration of the system is carried on directly in terms of y = In a where, 
we remind, a = a/a^e and a& e denotes the value of the scale factor when dark energy and 
matter give equal contribution to the total energy density. The variable y is directly related 
to the redshift so that by plotting the growth factor in terms of y also its redshift dependence 
can be easily sorted out since the following chain of equalities holds: 

y = In a = In (> de + 1) - \n(z + 1). (4.1) 

The physical range of y does not extend beyond 1, indeed the dependence of y max upon Wdc 
is monotonic and ranges between 0.338 (for Wdc = — 1 and for the parameters of Eq. (12.121) ) 
and 1.12 in the case w^e = —0.3 (where, strictly speaking the background geometry does not 
accelerate). It is practical to introduce the rescaled wavenumber 

„ 2 , k 2 ( k \ 2 (Zeq + 1) e (^c+i)y 



^ U 2 Uo ) ^MoOde + 1) [(*q + + Ode + l^-"^ + ^ cq + 1] ' ^ 

so that, for instance, Eq. (I2.43P will be rewritten as 

d5 c Idh N - d6 c - 

where, for a generic species X, 9 X = Ox/k. The evolution of £(£;, y) and h(k, y) is determined, 
respectively, by the following pair of equations 

+ [2 + i%)]^| = 2n R (R v <T u + R r cj B )-^{8 R + R i n B ) 
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[1 + w de )(c 2 de - w de )9 



dc 



d 2 h ,^ X1 dh 



6 fi R [<5 R + R^ B ] + 3 M 5 n 



dc 



;i + 34)5 de + 



9 



l + W de )(Cd e -W de )6' d e 



(4.4) 



(4.5) 



where Z(y) is expressible in terms of y by recalling the definition of Eq. (12. 6p and where the 
density contrasts in radiation and matter are defined as 



S R (k, y) = i? 7 5 7 (fc, y) + R v b v (k, y), 



S m {k,y) = —6 c {k,y) + —5 h (k,y). 



(4.6) 



Equations (14. 4 j) and (14. 5 p are derived as linear combinations of Eqs. (12.231) . ( 12. 36ft and ( 12. 3 7ft 
written in the y parametrization. Moreover, in Eqs. (14. 5p . (I4.4p and (14. 6 j) i? 7 = 1 — R u and 



Pu 3 x (7/8) x (4/11) 4 / 3 

p 7 + p, 1 + 3 x (7/8) x (4/II) 4 / 3 



0.4052 



(4.7) 



where 3 counts the massless neutrino families, (7/8) stems from the Fermi-Dirac statistics 
and (4/11) 4 / 3 comes from the kinetic temperature of neutrinos. Within the same notations, 
the evolution for the neutrinos (see Eqs. (I2.44p and (12.451) ) become: 



dy 
do v 
dy 



£(y) 



5 U - JC(y)a u , 



d5 u 
dy 



If ~ 

6 dy 3 



— + 6— 

dy dy. 



(4.8) 
(4.9) 



Prior to recombination the tight-coupling limit has been enforced. To zeroth order in the 
tight coupling expansion the photon quadrupole as well as the polarization vanish and the 
relevant evolution equations for the dipole and for the monopole will then become: 



de. 



7b 



+ 



Rb 



dy _Rb + 1 



-9 



7b 



£(y)&y 

A(R h + 1) A[R h (y) + 1] 



(4a B - Ob) 



d5 h ldh - d5 1 2dh 4 - 

^ = Ydy- KW ^ -dy- = 3dy--3 ,Ciy)e ^ 

After recombination the baryon and the photon velocities are different, i.e. 



(4.10) 
(4.11) 



(9b — 9~ 



— - ; zrzzr, , — fhH . ~ ; ; 



[l + R b (y)]]C(y)\ 



4Rb(y) 



+ 



, dO^ d9\) , 
1 dy dy ' 



(4.12) 



where the contribution of the quadrupole moment of the photon distribution has been ne- 
glected since it vanishes to zeroth-order in the tight coupling expansion; note that, as before, 
JC(y) = k/T-L while JC(y) = k/e'; c 2 h denotes the sound speed in the baryon-photon fluid 
J. The temperature of the baryons has been taken to coincide with the temperature 
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of the photons both before and after recombination. After a transient regime the evolution 
equations for the photons and for the baryons will obey separate evolution equations. More 
specifically the appropriate collective variables will be the total density contrasts of matter 
and radiation as well as the radiation and matter velocities: 

6 K {k,y)=R i e i {k,y) + R v e v {k,y), 9 m (k,y) = (— y) + (— W, v). (4-13) 

Finally the evolution equation for the dark energy (see Eqs. (12.471) and (I2.48P ) can be written 

as 

d5 dc 



dy 
+9< 



+ 3(c dc - w de )5 de + Ode + l)JC(y)6 dc 

1 + Wde dh 
2 ~d~y = 



1 + W de ) 2 - 

(c de - w dc )6 dc 



ic(v) 

dOde , n „ o \ 77 c 2 de IC{y) 

— h (1 - 3c de )V de + — —d dc = 0. 

dy (w dc + l) 



(4.14) 
(4.15) 



In the case of adiabatic initial conditions the asymptotic solution of the previous set of 
equations can be written as 



h(k,y) 



R 1 [4(7 B (k) - R„n B {k)) 4:R U + 5 



6(4R„ + 15) 



12(4R„ + 15) 



K.(k)\lC 2 (y), (4.16) 



K\y) 



ml - UR U - 75 
36 I '2(2.R I/ + 25)(4R I/ + 15) 



K*{k) 



+ 



i? 7 (15 - 20R P 



5 7 {k,y) 

$c(k,y) 
$b(k,y) 
5a e (k,y) 
6 yh (k,y) 

Uk,y) 
d c (k,y) 



10(4R„ + 15)(2.R I , + 25) 



[RM B (k)-4<r B (k)}jK, 4 (y), 
K 2 (y), 



-^— + a B {k) - — Q B {k) 



-R y n B (k) + - 



--RyQ B (k) + i 



n * (k) JC 2 (y), 



+ a B (k) - ^ B (k) 



~(i + w de )R,n B (k) + (Wdc 4 +1) ^(fc)/c 2 ( y ), 



-j-tt B (k) - a B 



-K*(k) + 



R u Q B {k) - Aa B {k) 



^a B (k) - ^Q B (k) 



£(v)- 



1 f (4/^ + 23) 



36 I 4R U + 15 



7e*(fc) 



= 0, 



(4.17) 
(4.18) 
(4.19) 
(4.20) 
(4.21) 
(4.22) 
/C%), (4.23) 

(4.24) 
(4.25) 
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Ode(k,y) 



dc 



4(2 -3c| e r 

-4w*o + 



4(2 - 3ci 
+ 



-K 3 (?/) 



(4.26) 



As a cross-check of the accuracy of the integration the Hamiltonian and the momentum 




. k= 0.003 Mpc -1 
(fully numerical) 

. k=0.003 Mpc" 1 
(analytical) 



Figure 4: In the left plot the growth factor is illustrated for w de = — 0.8 and in the case of 
three different wavenumbers, i.e. k = 0.3 Mpc -1 (full line), k = 0.03 Mpc -1 (dashed line) 
and k = 0.003 Mpc -1 (dot-dashed line). The right plot shows, for a particular wavenumber, 
the result of the numerical integration of this section and the analytical approximation based 
on Eq. fl3T32|) . 



constraints must always be satisfied and their explicit expression becomes, respectively, 



dh 

2/C(y)f - 3[fi R 5 R + £l M K + fide^de] = 0, 

ay 

d£ 3 3 

K(y)-T + 2Q R {y)6 R + -{w dc + i)Q de e dc + -o M i = o. 
ay 2 2 



(4.28) 
(4.29) 



The dimensionelss ratio defined in Eq. (14. 2 j) dominates against all the terms containing ex- 
plicitly the conductivity as implied by Eqs. ( I2.33p - (12.34j) . To get an idea of the inaccuracies 
introduced by neglecting the diffusivity terms and the terms containing powers of the electric 
field a further dimensionless ratio can be defined and it is given by K? a (y) = k 2 /(Ha). The 
ratio between JC a (y) and )C(y) gives 



£(v) 



1020.3 



(4.30) 



the inequality sign in Eq. (I4.30p arises since we assumed that in the dynamical range of 
the numerical integration f^Mo < Rm(z/)- Thus the inaccuracies introduced by the resistive 
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O -0.2 



B L = 5 nG 
B, = 10 nG 




Figure 5: The growth factor is illustrated in the for different values of la- 



terals are of higher order in comparison with accuracy goal and with the working precision 
of the numerical calculation. The variation of the parameters affecting directly the shape of 
the magnetized growth factor will now be illustrated. The wavenumbers directly probed by 
the currently available large-scale structure data range from /c min = 0.01 h Mpc -1 to k max = 
0.25/ioMpc -1 including also the typical range of scales at which the spectrum becomes 
nonlinear, i.e. k n \ > 0.09 h Mpc -1 . For illustration three different values of k are reported 
in the left plot of Fig. H i.e. k = 0.3 Mpc" 1 (full line), k = 0.03 Mpc" 1 (dashed line) and 
k = 0.003 Mpc -1 (dot-dashed line). Always in Fig. H] (plot at the right) the numerical 
calculation described in the present section is compared with the analytical results of section 
[3] for the same values of the parameters and by focusing the attention on the region of large 
y, i.e. y — > 1. For the parameters chosen in Figs. 0] and [5] the values of y at decoupling 
and for different barotropic indices are given, respectively, by ?/dcc = —6.569 for w = —0.8, 
Vdec = —6.616 for w = —0.9 and y^ec = —6.653 for w — — 1. The parameters chosen for the 
numerical integrations of Figs. H] and E] are the ones of Eq. ( 12. lip , where, however, w^e has 
been taken to be different from —1 to allow for dynamical dark energy fluctuations. The 
WMAP 7yr data alone are consistent with a rather broad interval of values of Wd e ranging 
from —0.63 to about —1.6. We shall not contemplate the possibility of having Wdc < —1: 
in the latter case sudden singularities can arise in the future. In Figs. HI [5] and [6] the 
speed of sound of dark energy is chosen to be c\ c = 1. The latter choice corresponds to 
identify the sound speed of dark energy with the speed of light as it happens, for instance, in 
quintessential models of dark energy. As already mentioned, the bounds on c\ e are not tightly 
fixed by observation and the differences induced by different values of c\ e will not be essential 
for the present discussion as long positive semidefinite. A comparison between the 

case = 1 and c^ e = is drawn in Fig. [6j In Fig. [5] the growth factor is illustrated for 
the variation of B^. The comparison of the numerical result with the analytical estimates 
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Figure 7: The variation of the sound speed of dark energy is illustrated for the two extreme 
values of the barotropic index adopted in the text. In the plot at the left Wde = —0.3; in the 
plot at the right w^ e = —0.9. 



is illustrated in the plot at the right. The differences between the numerical results and the 
analytical formula are determined, in this case, by the presence of dark energy fluctuations. 
Note that the range of y has been narrowed to emphasize the differences. The same aspect is 
illustrated in Fig. [6] where the growth factor is computed numerically in the case of different 
barotropic indices. In both plots of Fig. [6] the sound speed of the dark energy fluctuations 
coincides with the speed of light (i.e. c\ e = 1). The results presented so far support the view 
that general relativistic corrections alter (from 5 to 10 %) the growth rate from z ~ 100 
down to z ~ already for k ~ 0.003 Mpc -1 . 
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In the left plot of Fig. E]the value of the comoving wavenumber is fixed to k — 0.3 Mpc~ 
and the numerically computed gowth rates are illustrated for three different values of the 
barotropic indices of dark energy ranging between —0.3 up to —0.9. In the same figure (but 
in plot at the right) the numerical results are compared with the analytic approximations 
based on Eqs. ( I3.25P and f )3.32p -( 13.34p for a comoving wavenumber k = 0.03 Mpc -1 . A 
general feature of the obtained numerical results, as exemplified by Fig. E]is a disagreement 
between the analytic and the numerical results for very large or very small redshifts. For 
small y (i.e. large redshifts) the radiation must be properly treated (and it has been instead 
neglected in section [3]). Similarly, at low redshifts (i.e. large y) the fluctuations of the dark 
energy give a correction which has been neglected in the analytical discussion. In spite of 
the mentioned caveats, however, it seems that the analytical treatment correctly captures 
the region of intermediate redshifts and it is, overall, quite reasonable. While in the previous 
figures the sound speed of the dark energy has been taken to coincide with the speed of light, 
in Fig. [7J the variation of the sound speed is illustrated. In Fig. [7J in the left and right plots 
the barotropic indices take the two extreme values w dc = —0.3 and u>d e = —0.9. In each of 
the two plots two different values of the sound speed of the dark energy (i.e. c 2 dc = 1 and 
c de = 0) are illustrated. The range of redshifts has been narrowed to make the differences 
more visible. As expected the differences between the two cases are more evident in the 
region of low redshifts (i.e. large y), in agreement with the previous remarks. 

5 Concluding remarks 

Direct measurements of the growth rate of matter fluctuations can be used as a probe of 
large-scale magnetism. In the conventional ACDM paradigm (and in its extensions) the 
baryons fall into the dark matter potential wells so that their corresponding growth rate 
is determined primarily by the density contrast of cold dark matter. The latter statement 
must be partially revised in a magnetized environment where the resulting growth index 
is determined by the competition of the dark matter fluctuations and by the magnetic in- 
homogeneities entering the evolution equations of the baryons. The rich physical structure 
of the magnetized initial conditions of the Einstein-Boltzmann hierarchy has been explored 
in all phenomenologically interesting cases contemplating the standard adiabatic mode as 
well as the various isocurvature modes. The improved understanding of magnetized CMB 
anisotropies suggests a more thorough investigation of the effects of large-scale magnetization 
on the current paradigms of structure formation. The parameters minimizing the distortions 
of the temperature and polarization angular power spectra induce computable modifications 
on the shape of the growth rate. The magnetized growth rate has been scrutinized both for 
small and for large redshifts in comparison with the typical time-scale at which the baryons 
are freed from the drag of the photons. The results reported here pave the way for a more 
systematic scrutiny of the impact of large-scale magnetic fields on the current paradigm of 
structure formation based on the the ACDM scenario and on its neighboring extensions. 
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